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The past two decades of lunar exploration have seen the detection of 
substantial quantities of water onthe Moon’s surface. It has been proposed 
that a hydrated layer exists at depth in lunar soils, buffering a water cycle 
onthe Moon globally. However, a reservoir has yet to be identified for 

this hydrated layer. Here we report the abundance, hydrogen isotope 
composition and core-to-rim variations of water measured in impact 

glass beads extracted from lunar soils returned by the Chang’e-5 mission. 
The impact glass beads preserve hydration signatures and display 

water abundance profiles consistent with the inward diffusion of solar 
wind-derived water. Diffusion modelling estimates diffusion timescales of 
less than 15 years at a temperature of 360 K. Such short diffusion timescales 
suggest an efficient water recharge mechanism that could sustain the 

lunar surface water cycle. We estimate that the amount of water hosted by 
impact glass beads in lunar soils may reach up to 2.7 x 10“ kg. Our direct 
measurements of this surface reservoir of lunar water show that impact 
glass beads can store substantial quantities of solar wind-derived water on 
the Moon and suggest that impact glass may be water reservoirs on other 


airless bodies. 


It has long been argued that there could be water and other volatile 
species at the surface of the Moon’ °. Renewed lunar exploration and 
advances in remote-sensing measurements in the late 1990s allowed 
the neutron spectrometer on board the Lunar Prospector mission to 
confirm the existence of water ice at the lunar poles*. Following this, 
the Moon mineralogy mapper instrument on board the Chandrayaan-1 
spacecraft detected the diagnostic absorption bands of hydroxyl 


and/or water at 2.8-3.0 um onthesurface of the Moon’. Detailed analyses 
indicated that hydroxyl/water could be present all over the Moon’s 
surface, which also displayed temporal and spatial variations, with 
equivalent water abundances ranging between -10 and 1,000 ug g“ 
(refs. €). Furthermore, the Lunar Crater Observation and Sensing 
Satellite impact experiment carried out in 2009 provided direct evi- 
dence for high water-ice abundances (5.6 + 2.9 wt%) in permanently 
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Fig. 1| Representative back-scattered electron images of CE5 impact 

glass beads after nanoSIMS analyses. a, CE5#33,002. b, CE5#33,003. 

c, CE5#33,036. d, CE5#33,046. e, CE5#33,052. f, CE5#33,076. These impact glass 
beads have spherical shapes and have homogeneous chemical compositions 


S28 of bees BS 


Profile 2 


g aljold 


vu 
ss 
° 
=a 
© 
= 


Profile 5 l 


(Supplementary Table 2). Beads CE5#33,003 and CE5#33,036 are coated by finer 
agglutinates. Six profiles carried out on five impact glass beads in the second 
session are outlined in red. Blue squares indicate locations of NanoSIMS in the 
first analytical session. Each NanoSIMS analysis pit is 7 x 7 um’. 


shadowed regions within Cabeus crater". Elevated water-ice abun- 
dance in lunar polar regions was further supported by the neutron 
flux measurements performed by the Lunar Exploration Neutron 
Detector on board the Lunar Reconnaissance Orbiter spacecraft”. 
Recently, the neutral mass spectrometer on the Lunar Atmosphere 
and Dust Environment Explorer detected exospheric water liberated 
by meteoroid impacts”, and ground-based telescope observations 
detected molecular water on the lunar surface". Today, there is little 
doubt that most of the Moon’s surface harbours water in one form or 
another. However, the origin(s) of this lunar surface water and its spatial 
distribution and evolution during regolith gardening remain largely 
unknown, despite key implications for future lunar surface exploration 
and for better understanding the (sub)surface water reservoir and 
processing on Solar System airless bodies. 

There are several potential sources and processes that could 
have contributed to the water inventory at the surface of the Moon”, 
suchas (1) solar wind implantation, (2) outgassing of volatiles dur- 
ing lunar volcanism” ~’, (3) deposition of volatile-bearing pyroclastic 
deposits” and minerals” and (4) delivery by impacts of comets and 
asteroids”*”’. It is generally thought that solar wind hydrogen-ion 
implantation could react with surface minerals to produce hydroxyl or 
water in lunar soils*°’. The surface water produced and/or delivered 
onthe Moon’s equatorial regions may migrate to polar regions, driven 
by temperature oscillations”, and could also be partially released 
to space”. Hence, a scenario of lunar water cycle was further proposed 
to describe the retention, release and replenishment of water on the 


surface of the Moon”. To sustain a water cycle at the surface of the 
Moon, there should be a hydrated layer (reservoir) at depth in lunar 
soils’. However, finding this water reservoir has remained elusive, 
despite several studies having investigated the water inventory of fine 
mineral grains in lunar soils (for example, refs. ***°), impact-produced 
agglutinates”, volcanic rocks (for example, ref. *°) and pyroclastic 
glass beads”"*°. Consequently, there must be a yet-unidentified water 
reservoir in lunar soils that has the capacity to buffer a lunar surface 
water cycle. 

The lunar soils returned by the Apollo, Luna and Chang’e-5 (CE5) 
missions comprise mainly rock and mineral fragments, pyroclastic 
glass beads, impact-produced agglutinates and impact glass beads”. 
Among these components, the water inventory of impact glass beads 
has not yet been investigated in detail, despite these glassy beads being 
potential candidates for playing a significant role in a lunar surface 
water cycle. To investigate this possibility, we carried out asystematic 
characterization of the petrography, major element composition, 
Raman characteristics, water abundance and hydrogen isotope com- 
position on the impact glass beads returned by the CES mission, aiming 
to identify and characterize the missing water reservoir on the Moon’s 
surface. A recent geochronological study of CE5 impact glass beads has 
shown that they formed more or less continuously for the past 2 Gyr, 
with prominent peaks in formation ages at ~575 million years ago (Ma), 
380 Ma, 68 Ma and 35 Ma, originating predominantly from a few 1- to 
5-km-diameter impact craters in the Em4 (Eratosthenian-aged mare) 
basaltic unit where CE5 landed”. 


Nature Geoscience 


Article 


https://doi.org/10.1038/s41561-023-01159-6 


Water abundance and hydrogen isotope 
composition of impact glasses 

The studied impact glass beads were hand picked from a lunar 
soil sample (sample CESCO1O00YJFMO0103; ~1 g) scooped by the CE5 
lander robotic arm. A total of 117 individual spherical glass beads were 
characterized using field-emission scanning electron microscopy 
(SEM), electron probe microanalyser (EPMA) and Raman spectro- 
scopy (Fig. 1, Supplementary Figs. 1 and 2, Supplementary Table 1 
and Extended Data Fig. 1). Among them, 32 impact glass beads charac- 
terized by smooth exposed surfaces and chemical compositions 
consistent with those of bulk CES basalts were chosen for in situ water 
abundance and H isotope analysis using a nanoscale secondary ion 
mass spectrometer (NanoSIMS) (Fig. 1, Extended Data Figs. 2 and 3 
and Supplementary Table 1). The CE5 impact glass beads con- 
tain 0-1,909 ug g” equivalent H,O, with H isotope compositions 
(given in delta notation, 6D = 1000 x (D/Hgampie/D/Hsmow-1), where 
D/Hgampie is the measured D/H ratio measured in the sample and D/ 
Hemow is the D/H ratio (1.5576 x 10“) of the standard mean ocean water) 
ranging from —990 + 6%o to 522 + 440%o and which are negatively 
correlated with the water abundances (Fig. 2). We also analysed water 
abundances and H isotope compositions along six transects in five 
glass beads (CE5#33,003, CE5#33,036, CE5#33,046, CE5#33,052 and 
CE5#33,076; Fig. 1), which show elevated water abundances at the rims 
that gradually decrease towards the cores while 65D values decrease 
from cores to rims (Fig. 3 and Extended Data Figs. 4 and 5). Three glass 
beads display sharp decreases in water content from the rims towards 
the cores, while the other two glass beads display a smoother water 
abundance variation (Fig. 3 and Extended Data Figs. 4 and 5). 


Solar wind origin of water 

The rims of CES impact glass beads are generally characterized by 
higher water abundances (up to -2,000 ug g”) and lower SD values 
(--990%o) compared with pyroclastic glasses, impact-produced agglu- 
tinates and melt inclusions reported in previous studies (Fig. 2 and 
Supplementary Table 4). The extremely low &D values measured at 
the rims of CE5 impact glasses (Fig. 3) are consistent with the H iso- 
tope composition of solar wind*°“". These low 8D values are distinct 
from other potential lunar water reservoirs, such as indigenous water 
outgassed during lunar volcanism (6D = -300 to 1,200%v) (refs. °°"), 
water-rich carbonaceous chondrites (65D = -200 to 800%v) (refs. ****) 
and comets (6D =0 to 2,500%v) (refs. ***°), indicating that the water 
retained in CE5 impact glasses is of a solar wind origin. 

There are two potential mechanisms, de-gassing and two- 
endmember mixing, that could explain the negative correlation 
between water abundances and SD values defined by CES impact 
glasses (Fig. 2). Because the initial D/H ratio is so low (D/H =3.1 x 10° 
for the analysis with the highest H,O abundance), the H, de-gassing 
modelling indicates that D/H fractionation due to H, loss would yield 
an increase in 6D values from about -980% up to only -926%o for 99.9% 
loss, starting with an initial water abundance in CE5 impact glasses 
of -2,000 ug g” (Fig. 2 and Methods). Therefore, H, loss from 
CE5 impact glasses via de-gassing is not a viable mechanism to 
account for the observed SD variations from around -990%o to 
+522%o (Fig. 2). More important, the high water content at the rim 
but low water content in the core of the glass beads is not consistent 
with the de-gassing scenario. 

In the other scenario, the CES impact glasses could have lost 
most of the original water inventory of their precursors during their 
formation at high temperature, following which solar wind-derived 
water started diffusing inwards into these glass beads on the Moon’s 
surface. This scenario is consistent with the notable lower abun- 
dances of volatile elements such as Na,O (most <0.05 wt%) and K,O 
(<0.06 wt%) in CES impact glasses compared with CES mare basalt 
fragments (Na,O = 0.6 wt% and K,O = 0.1-0.2 wt%, on average)”. This 
binary mixing model, between a water-poor but relatively high D/H 
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Fig. 2 | Water abundances and ôD values measured in CE5 impact glass 

beads. All glass bead analyses define a negative correlation between the water 
abundances and SD values, which is consistent with a two-endmember mixing 
model (red lines) that indicates addition of solar wind-derived water with 6D 
value of -990%o into a low water abundance (H,0, = 5-50 pg g”) CE5 impact glass 
with 8D value of +500%o (ref. °°). De-gassing of H, from CE5 impact glasses cannot 
be responsible for the large observed H,O-SD variations (brown line). The inset 
diagram is the zoom-in view of water abundance range of 0-100 ug g”. For more 
details, see Methods. The literature data for melt inclusions, pyroclastic glasses 
and impact-produced agglutinates are also plotted for comparison, and the 
datasets are listed in Supplementary Table 4. All error bars correspond to 2 s.d. 
analytical uncertainties. SMOW, standard mean ocean water. H,O, is the water 
abundance of impact glass beads prior to the inward diffusion of solar wind- 
derived water. 


endmember and a water-rich but extremely low D/H endmember, 
reproduces well the range of water abundances and 8D values meas- 
ured in CE5 impact glasses (Fig. 2). Modelling requires that the initial 
water abundances of CE5 impact glass beads are less than 50 pg g” 
before the addition of solar wind-derived water (Fig. 2). This require- 
ment is consistent with the water abundances measured at the cores 
of CE5 impact glasses, which are comparable to the instrument H,O 
background (-10-30 ug g”; Methods), indicating that any water poten- 
tially present in the precursor materials was lost during formation of 
the CE5 impact glass beads, before the addition of solar wind-derived 
water. Most important, the hydration profiles in CES impact glass 
beads show that solar wind-derived water has diffused inwards 
into the glasses (Fig. 3), unambiguously demonstrating that the 
addition of solar wind-derived water post-dated the formation of the 
impact glass beads. 


Impact glasses as the driver for the lunar surface 
water cycle 

The negative correlation between water abundances and SD values, 
coupled with the hydration profiles observed for CE5 impact glass 
beads, reveal that solar wind-derived water can diffuse into and be 
stored in lunar impact glasses on the Moon’s surface, probably via 
post-implantation diffusion’’”’”. To model the diffusion timescale, we 
used Raman spectroscopy to confirm that water hosted in CE5 impact 
glass beads is inthe form of hydroxyl/molecular water (Extended Data 
Fig. 6). The diffusion coefficient of hydroxyl/water in glasses** was 
used for modelling the duration timescale necessary for forming the 
observed hydration profiles of CE5 impact glass beads (Methods). 
Diffusion modelling results indicate that the diffusion time needed 
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Fig. 3 | Two profiles across CE5 impact glass bead CE5#33,036. a,b, Water 
abundance (a) and 8D values (b) measured along Profile 1, as outlined in Fig. 1c. 
c,d, Water abundance (c) and 5D values (d) measured along Profile 2, as outlined 
in Fig. 1c. The apparent distances along the profiles were measured from the rim 
on aside (Fig. 1c). A water-diffusion coefficient of 20.84 um? yr ‘at T= 360 K was 
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used**. The maximum diffusion durations for Profile 1 and Profile 2 are 9 and 15 yr 
at T=360 K, respectively. Profile 2 was further modelled by a post-diffusion 
de-gassing on the basis of the bell shape of water abundance on the right side. 
More modelling details can be found in Methods. All error bars correspond to 
2s.d. analytical uncertainties. 


to produce the observed hydration profiles varies from 1 year to 
15 years at a peak temperature (7) of 360 K, relevant forthe landing site of 
CE5” (Fig. 3, Extended Data Figs. 4 and 5 and Methods). This short 
diffusion time indicates that the solar wind-derived water can be 
rapidly accumulated and stored in lunar impact glass beads. It is inter- 
esting to note that one profile on bead CE5#33,036 displays notable 
post-diffusion water loss at the very rim (Fig. 3, and Methods), suggest- 
ing that CES impact glass beads have the capacity of both incorporat- 
ing solar wind-derived water and releasing it to the lunar exosphere 
if T increases. Therefore, we propose that impact glass beads in 
lunar soils are a prime water reservoir candidate able to drive the 
lunar surface water cycle (Fig. 4). 

Linking a lunar surface water cycle with impact glass beads 
implies three major processes: (1) formation of impact glass beads, 
(2) solar wind-derived water implantation and diffusion into the 
beads and (3) water release to the lunar exosphere during Tincrease, 
possibly due to sunlight irradiation and/or meteoroid impacts 
(Fig. 4). Impact glass beads are produced during melting of lunar 
soil and regolith targets when impactors hit the Moon. Most of the 
water present in the precursor materials would have probably been 
lost during high-temperature melting that led to the formation 
of the impact glass beads. Solar wind hydrogen ions would then 
be continually implanted into lunar soils and impact glass beads 
and would combine with oxygen atoms to form the hydroxyl/water 
that diffused into the glass interiors and is stored in glass beads. This 
water stored in impact glass beads could then be released to the lunar 
exosphere due to sunlight irradiation” and/or meteoroid impacts”. 


Some of the CE5 impact glass rims contain higher water abun- 
dance (up to -2,000 ug g`”) than that reported in Apollo agglutinates 
(<600 ug g’ (Fig. 2); however, they may still not be water saturated 
compared with terrestrial basaltic glasses of similar chemical compo- 
sitions that can contain up to 2.43 wt% H,O (ref. °). This could be the 
result of dynamic diffusion and release of water in the impact glass 
beads controlled by the time-of-the-day temperature oscillations”. 
The dynamic ingress and egress of water in impact glass beads could 
have acted as a buffer to explain the global and daily variations of water 
abundance on the lunar surface and in the lunar exosphere”*”. 


Inventory of water in lunar soils 

Using the water abundance measured in CES impact glass beads and 
mare basalts”, we can estimate the relative contribution of solar wind 
implanted in impact glasses and basalts to the global water inventory in 
the CES regolith. Considering the zoning features of water in CES impact 
glass beads, we estimate that these beads have a bulk water abundance 
of 132-1,570 ug g”, which is notably higher than that of the bulk CE5 
mare basalts (7 + 3 pg g”) (ref. °°). The modal abundance of impact glass 
beads in lunar soils is inthe region of 3-5 vol.% estimated from Apollo 
lunar soils”. Using this modal abundance yields a contribution by the 
CE5 impact glass beads to bulk soils of 4-78 ug g” equivalent water, 
whichis consistent with bulk water abundance estimates derived from 
reflectance spectroscopy measurements carried out from the CES 
lander”. Importantly, the major water-bearing phase in CE5 basalts is 
apatite in which water occurs as hydroxyl occupying the X site in apatite 
crystal structure****. The daily temperature variation of 93 K to 423 K 
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Fig. 4 | Schematic diagram of the lunar surface water cycle associated with 
impact glass beads. a, Because of the high temperature of formation ofimpact 
glass beads, most water present in precursor materials would have been lost. 

b, After their deposition at the surface, solar wind-derived water would diffuse 
into impact glass beads after solar wind H* implantation. c, Gardening would 
transfer the impact glass beads deeper into the soils, creating a water reservoir 
at the subsurface of the Moon. The impact glass beads at the surface would still 
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be able to release water into the lunar exosphere, due to meteoroid impact, 
for example”. The amount of water held by impact glass beads in lunar soils 
is estimated to be 0.03-27.00 x 10” kg, assuming a global depth of lunar soils 
of 3-12 m (refs. °°’), a modal abundance of impact glass beads in lunar soils 
of 3-5 vol.% (ref. °') and impact glass beads having a bulk H,O abundance of 
132-1,570 pg g”. 


onthe Moon’s surface is not high enough to drive hydroxyl out from 
the apatite structure’, indicating that impact glass beads are 
probably the dominant reservoir involved in the lunar surface water 
cycle, except for the space-weathered rims and possibly water ice in 
polar regions. 

The global lunar surface is covered by lunar soils with a thickness 
ranging from about3 mto12 m (refs. °°’). The amount of water in lunar 
soils contributed by impact glass beads varies from 3.0 x 10" kg to 
2.7 x 10“ kg, based on the bulk water abundance of 132-1,570 pg g'in 
impact glass beads and a modal abundance of -3-5 vol.% (ref. ") (Fig. 4). 
Importantly, this estimate provides a minimum quantity of water avail- 
able at the lunar subsurface as it does not take into account the water 
contribution from permanently shadowed regions* ©", indigenous 
lunar rocks (for example, ref. 7°), delivery by meteoritic impacts?” 
and solar wind-derived water preserved in the space-weathered rims 
of lunar soils****. These findings indicate that the lunar soils contain 
a much higher amount of solar wind-derived water than previously 
thought, which could be a water reservoir for in situ utilization in 
future lunar exploration. Indeed, this water entrapped in impact glass 
beads appears to be quite easy to extract as shown by the post- 
diffusion de-gassing profile for bead CE5#33,036 (Fig. 3 and Methods). 


Implications for water on airless bodies 

Asteroid and comet impacts are the major exogenous processes that 
reshape the surface morphologies of airless bodies, as evidenced by the 
widespread presence of impact craters on the Moon”, Mercury” and 
asteroids“. These impacts probably create impact glasses and glass 
beads on any airless bodies, as evidenced by the studies of lunar soils 
and howardites®. This study demonstrates that impact glass beads have 
the capacity to store significant quantities of solar wind-derived water 
at the surface of airless bodies, in addition to the possible presence of 
water ice trapped in permanently shadowed areas in polar regions. The 
presence of water, stored in impact glass beads, is consistent with the 
remote detection of water at lower-latitude regions of the Moon’ ”, 
4 Vesta™ and Mercury°***. Our findings indicate that the impact 
glasses on the surface of Solar System airless bodies are capable of 
storing solar wind-derived water and releasing it to space. 


Online content 

Any methods, additional references, Nature Portfolio reporting sum- 
maries, source data, extended data, supplementary information, 
acknowledgements, peer review information; details of author con- 
tributions and competing interests; and statements of data and code 
availability are available at https://doi.org/10.1038/s41561-023-01159-6. 
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Methods 

Sample preparation 

The CE5 lunar soil allocated by the China National Space Administration 
and used inthis study is CESCO100YJFMO0103, weighing ~1 g. The CES 
lunar soils were scooped by the robotic arm on board the CES lander. 
A total of 150 individual grains, varying in size from ~50 um to -1 mm 
in diameter, were hand picked under a binocular microscope in the 
ultraclean room at the Institute of Geology and Geophysics, Chinese 
Academy of Sciences (IGGCAS). The hand-picked grains were mounted 
in epoxy and prepared as a double polished section witha thickness of 
about 100 um (Supplementary Fig. 1and Supplementary Table 1). The 
surrounding epoxy was removed to improve the vacuum quality for 
the water abundance measurements using NanoSIMS”. The prepared 
section was cleaned using anhydrous ethanol before drying at 60 °C 
in a baking oven overnight. 


SEM observation 

Petrographic observations were carried out using field-emission SEM 
using FEI Nova NanoSEM 450 and Thermofisher Apreo instruments 
at the IGGCAS, using electron beam currents of 2.0-6.4 nA and an 
acceleration voltage of 15 kV. The prepared sections were initially 
coated with carbon for petrographic observations and analyses of 
mineral chemistry and then recoated with Au for in situ NanoSIMS 
measurement of water content and hydrogen isotopes. After NanoSIMS 
analyses, Au coating was removed, and the samples were recoated 
with carbon and observed by SEM to confirm the positions of the 
NanoSIMS spots. 


Electron probe microanalysis 

Before the NanoSIMS analyses, we used the CAMECA SXFive EPMA 
at IGGCAS to quantify the major and minor elemental abundances 
in CES impact glass beads. The operating accelerating voltage was 
20 kV and the beam current was 10 nA, with a 10 um beam diameter. 
The total counting time was 10 min for each analysis. One basalt glass 
(VG-2) and one komatiite glass (MPI-DING-GOR128) were used to 
monitor the analysis accuracy and precision®. The detection limits 
were (30) 0.02 wt% for Na, Ni, K, S and P, 0.04 wt% for Ca, 0.05 wt% 
for Al, Mg and Cr, 0.06 wt% for Si and Ti and 0.09 wt% for Mn and 
Fe. The EPMA data obtained for impact glass beads are listed in 
Supplementary Table 2. 


In situ water abundance and hydrogen isotope analyses 

Two analytical sessions were designed to measure the water abun- 
dances and hydrogen isotope compositions of CE5 impact glass beads 
using a CAMECA NanoSIMS SOL at the IGGCAS. In the first session, 
measurements at the cores and rims of CE5 impact glass beads were 
carried out to investigate the potential water zonation features in 
the glasses. In the second session, more-detailed profile measure- 
ments on five individual CE5 impact glasses were acquired, as well as 
several parallel measurements on other glass beads for cross checking. 
A liquid nitrogen cold trap was used to further improve the vacuum 
quality in the second session. 

The sample and the standard prepared in different mounts were 
loaded in the same sample holder of NanoSIMS and were baked over- 
night at ~60 °C inthe NanoSIMS airlock. The holders were then stored 
in the NanoSIMS sample chamber to improve the vacuum quality and 
minimize the H background**”. The vacuum pressure in the anal- 
ysis chamber was 1.5 x 10™ to 1.9 x 10°” Torr during analysis. Each 
7 um x 7 um analysis area was pre-sputtered for 2 min with a Cs* ion 
beam current of 2 nA to remove the surface coating and potential 
contamination. During analysis, the secondary anions !H_, 2D", 2C~ 
and O~ were simultaneously counted by electron multipliers (EMs) 
from the central 5 um x 5 um areas using the NanoSIMS electronic 
gate technique (49% blanking). A 44 ns dead time was corrected for all 
EMs, while the EM noise (<107 cps) was ignored. We used a primary ion 


beam current of -0.5 nA for analysis, corresponding to a beam size of 
~500 nm in diameter. The charging effect on the sample surfaces was 
compensated by an electron gun during analysis. 

A chip of the nominally anhydrous San Carlos olivine reference, 
with a reported water content of 1.4 ug g” (ref. °°), was used for 
instrument H,O background corrections, following the relationship 
H/Opg = (Heounts ~ Hg)/Ocounts and D/Hyreasured = (1 —f) X D/Hirue +f D/H, 
where fis the proportion of H emitted from the instrumental back- 
ground”. In the first session, D/H,, was 2.83 (+1.98) x 10* and 
Hyg = 815 + 457 cps (2 s.d., N=16, corresponding to an instrument 
H,O background abundance of 29.8 + 14.8 ug g” (2 s.d.)). In the second 
session, the measured H counts on the central areas of CE5 impact 
glass beads are even lower than that of San Carlos olivine. Therefore, 
the extremely low H/O ratio measurements on CES impact glass 
beads were treated as reference (H,, = 215 + 40 cps, corresponding 
to an instrument H,O background of 11.7 + 2 ug g”) for background 
correction. 

After an instrument H,O background subtraction, the water 
abundances of CE5 impact glass beads were calculated from the 
background-subtracted H/O ratios multiplied by the slope of the cali- 
bration lines (Extended Data Fig. 2), which were determined by meas- 
uring two apatite standards, Durango apatite (H,O = 0.0478 wt% and 
&D = -120 + 5%o) (refs. 7°”) and Kovdor apatite (H,O = 0.98 + 0.07 wt% 
and &D = -66 + 21%) (ref. ”), the SWIFT mid-ocean-ridge basalt (MORB) 
glass (H,0 = 0.258 wt% and 6D = -73 + 2%o) and two basaltic glasses, 
519-4-1 (H,O = 0.17 wt%) (ref. °) and 1833-11 (H,O = 1.2 wt%) (ref. *’) 
(Supplementary Table 3). Corrections for instrumental mass frac- 
tionation (IMF) on H isotopic compositions were carried out using 
repeated analyses of the SWIFT MORB glass standard and monitored 
by analysing both the Durango apatite and Kovdor apatite standards 
during the whole analytical session (Extended Data Fig. 3). The matrix 
effects on water abundance and IMF on hydrogen isotope composi- 
tion are the same between apatite and silicate glass within analytical 
uncertainties”. Hydrogen isotopic compositions are given using the 
delta notation, 8D = ((D/H) sampie/ (D/H) smow) ~ 1) x 1,000%, where SMOW 
is the standard mean ocean water with a D/H ratio of 1.5576 x 10. 
More technical details can be found in ref. ”. All data are reported with 
their 2 s.d. analytical uncertainties, which include reproducibility of 
D/H measurements on the reference materials, uncertainty of 
H,O background subtraction and internal precision on each analysis 
(Supplementary Tables 3 and 4). The raw measured D/H ratios were 
corrected for the background, followed by correction for IMF. 


Raman spectroscopy 

The Raman spectra of impact glasses were collected using a WITec 
alpha300 R confocal Raman system and carried out after the NanoSIMS 
measurements at the IGGCAS. An optical microscope with an objective 
of 0.9 NA and x50 magnification was used to focus the Ar’ laser beam 
with a 484 nm excitation wavelength on the target phases. Raman 
spectra were acquired with a total integration time of 500 s, a 9.2 mW 
laser power and a beam size of 5 um. Three CES impact glass beads 
CE5#33,076 with water abundances zoning from1,909 ug gat the rims 
to <10 ug g ‘in the cores were chosen for Raman analyses (Extended 
Data Fig. 6). Along with the CE5 impact glass beads, terrestrial MORB 
glasses, MRN-G1 (H,O = 1.5 wt%) (ref. ”) and EPR-G3 (H20 = 0.22 wt%)”, 
were also measured to demonstrate the variation of Raman peak inten- 
sity of water and/or hydroxyl at -3,300 cm“ to 3,700 cm” with their 
water abundance (Extended Data Fig. 6). 


Correction of water abundances and D/H ratios for spallation 
effects 

The measured D/H ratios have been corrected for the potential 
effects of spallation by cosmic ray, using a D production rate of 
2.17 x 10” mol g” Myr” (ref. *) for impact glasses. The correction errors 
induced by D spallation are around 5% on D/H ratios and negligible 
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on water content”. Because no cosmic-ray exposure age is yet avail- 


able for the CES impact glass beads, a cosmic-ray exposure age of 
50 Ma estimated by Hu et al.°° was used for correction of spallation 
in this study (Supplementary Table 3). Another D production rate, of 
9.2 x 10° mol g? Myr” (ref. ”), was also used for correction of spalla- 
tion effects to account for the suggestion that the D production rate 
of 2.17 x 10°? mol g” Myr” (ref. “) may not be appropriate for impact 
glasses. As the D production rate in Ref. ” is slower, corrections for 
spallation-produced D are even smaller, not making any notable dif- 
ference overall for our analytical results (Supplementary Table 4). 


De-gassing modelling 

The hydrogen isotope fractionation during volatile loss into a 
vacuum is given by a” = M1/M2, where M1 and M2 are the masses of 
the volatile phase isotopologues. The change of the isotopic 
composition of H during volatile loss by Rayleigh fractionation is 
given by R = R, xf“ », whereR, and Rare the initial and final D/H ratios 
for a fraction fof remaining hydrogen”. De-gassing of H, (M1 =2 for 
H, and M2 =3 for HD) yields an a value of -0.8165 (ref. ”°) (Fig. 2). 
The de-gassing modelling assumes that the CES impact glass beads 
have an initial water abundance of 2,000 ug g”, an initial D/H ratio 
of 3.1 x 10% (corresponding to that of the analysis with the highest 
H,O abundance) and de-gas in the form of H,, yielding a 6D value 
increase from -980% at f= 1 to —926%o at f= 0.001. 


Diffusion modelling 

The five sharp hydration profiles of CES impact glass beads were 
simulated by a simplified diffusion process that follows Fick’s second 
law as equation (1) (ref. ””). 


x 
2VDt 


C(x) = (Cs — CO) erfe( ) + CO (1) 


In equation 1, D is the diffusion coefficient that depends on tempera- 
ture, xis the distance to the rim, CO is the initial water abundance of 
the CE5 impact glass beads, C(x) is the water abundance at distance x, 
Cs is the maximum water abundance at the rims, t is duration of 
diffusion and erfc is the complementary error function. 

Considering the scooped samples were collected from less than 
3 cm deep”, diffusion simulations were performed using a peak tem- 
perature of 360 K, relevant fora latitude of 45° N on the Moon”, because 
the peak temperature is the main controlling factor for the diffusion 
time”. The H diffusion coefficients of water/hydroxyl (D,,) determined 
on silica glass*’, D,, = 20.84 jm? yr‘ at 360 K, were used in our simula- 
tions. It is not possible to reconstruct the orientation of the exposed 
surface of the studied impact glass beads on the lunar surface. Thus, 
we carried out a simplistic diffusion modelling based on the measured 
profiles, assuming a constant Cs at the rims for each bead based on their 
diverse Pb-Pb ages”. The diffusion duration is roughly a maximum 
constraint for qualitatively estimating the efficiency of solar wind- 
derived OH/H,O recharged in impact glass beads. The water abun- 
dances at the rims (Cs) of the six CE5 impact glass bead profiles 
are extrapolated from the measured profiles using 2,000 pg g7, 
2,000 pg g”, 1,600 ug g“, 500 pg g“, 2,600 ug g” and 3,000 pg g7 
and a CO of 10 pg g` for Profiles 1-6, respectively. The diffusional 
durations vary from 3.0 to 9.0 years, 5.0 to 15.0 years, 1.0 to 4.5 years, 
2.0 to 6.0 years, 0.5 to 4.0 years and 2.0 to 14.0 years at T= 360K 
for Profiles 1-6, respectively, based on the apparent distance (Fig. 1) 
and a fixed Cs at the rims for each profile. 

For the notable post-diffusion de-gassing feature of Profile 2, a 
two-stage modelling of water diffusion (Stage 1) and post-diffusion 
de-gassing (Stage 2) was carried out to fit the measured water abun- 
dance profile. The water-diffusion details in Stage 1 have been described 
in the preceding. A post-diffusion de-gassing of water (Stage 2) was 
modelled to fit the bell shape of water abundance on the right side 


(Fig. 3). The water abundance on the right side modelled in Stage 1 was 
taken as the initial condition for modelling of Stage 2. The de-gassing 
rate is assumed to be proportional to the surface water abundance 
along the profile, where transmission capacity of a boundary was 
formulated as an additional boundary condition for the diffusion 
equation (see ref. * for more details). The corresponding 6D profiles 
are calculated by a binary mixing model (Two-endmember mixing 
modelling) considering a limited diffusional fractionation for H iso- 
tope” compared with the large difference in 6D values for the two 
endmembers (8D = —980%o, D/H = 3.0 x 10° for solar wind-derived 
water and 5D = 500%o, D/H = 230 x 10° for the initial water in 
impact glass beads). 


Two-endmember mixing modelling 

In the two-endmember mixing model, the solar wind endmember 
is set to have a 6D value of —980%o (ref. ““) and a water abundance of 
2,000 ug g”, which corresponds to the maximum measured value 
in this work (Fig. 2). The initial water component in CE5 impact glass 
beads could havea minimum SD value of ~500%v based on the measure- 
ments on CE5 apatite”, yielding a best fit of initial water abundance 
of 5-50 ug g` before the addition of solar wind-derived water (Fig. 2). 
On the basis of the following mass balance (equations (2) and (3)), we 
can model the relationship between the measured water abundances 
and ôD values: 


H200 + H2O solar = H2 Omeasured (2) 


H200 x D/Ho of H2Osolar x D/Hsotar = H2Omeasured x D/H measured (3) 


where HO, and D/H, are the initial water abundance and D/H ratio 
of CE5 impact glass beads before the addition of solar wind-derived 
water, respectively, H,O,,),, is the water abundance contributed by 
solar wind-derived water, D/H,1., is the hydrogen isotopic composi- 
tion of solar wind-derived water, and H30 measurea ANd D/H measured are 
the water abundance and D/H ratio measured from CE5 impact glass 
beads, respectively. 


Petrography and chemistry of CE5 impact glass beads 
Petrography. A total of 150 impact glass fragments were mounted in 
CE5#33 mount (Supplementary Fig. 1 and Supplementary Table 1). 
There are 117 glass fragments displaying ovoid to spherical shapes, 
which are mainly composed mainly of glass, 1 glass rod, 3 olivine and 
plagioclase monomict clasts, 1 basaltic fragment and 28 glass frag- 
ments showing irregular shapes, which are composed mainly of glass, 
mafic mineral fragments and large bubbles (Supplementary Fig. 1 
and Supplementary Table 1). The ovoid to spherical fragments are 
referred to as glass beads in this work and further classified into two 
textural types, homogeneous and heterogeneous glass beads (Fig. 1, 
Supplementary Fig. 2and Supplementary Table 1). The homogeneous 
glass beads display smooth surfaces and homogeneous contrast in 
back-scattered electron images while the heterogeneous glass beads 
usually contain some vesicles, pyroxene and plagioclase clasts and 
micrometre-sized sulfide spheres (Fig. 1 and Supplementary Fig. 2). 
The textural differences between homogeneous and heterogeneous 
glass beads could reflect a higher formation temperature for the 
homogeneous glass beads (Fig. 1 and Supplementary Fig. 2). 


Chemistry. A total of 32 homogeneous impact glass beads were 
selected for EPMA analyses. The analytical results are listed in 
Supplementary Table 2. The analysed glass beads have a mare 
origin, except for two glass beads plotting in the highlands popula- 
tion and another two distributed in pyroclastic areas judging from 
their CaO/AI,O, and MgO/AI,O, ratios (Extended Data Fig. 1)’”°°. These 
glasses plotting in non-mare fields are not discussed further in this 
work because their origin is currently uncertain. The homogeneous 
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CE5 impact glass beads are composed mainly of SiO, (31.2-44.5 wt%), 
TiO, (3.95-6.54 wt%), Al,O; (10.4-17.3 wt%), FeO (17.6-23.5 wt%), 
MgO (5.94-10.9 wt%) and CaO (10.1-15.4 wt%), with minor Cr,0, 
(0.12-0.28 wt%), MnO (0.22-0.31 wt%), NiO (0.15-0.61 wt%), Na,O 
(<0.17 wt%, most <0.05 wt%) and K,O (<0.06 wt%) (Supplementary 
Table 2). Each single glass bead is homogeneous in terms of its major 
element chemical composition (Supplementary Table 2). Most CES 
homogeneous impact glass beads have compositions similar to the 
bulk CE5 basalt compositions reported by refs. ***", indicating that they 
could have been formed locally rather than ejected from other regions. 
This prediction is consistent with the similar geochemical characteris- 
tics between CE5 impact glass beads and CE5 bulk basalts”. The Pb-Pb 
dating of the CES impact glass beads indicates that they were formed 
in the range of 4-2,000 Ma by various local meteoritic impacts”. 


Water abundances and hydrogen isotope compositions of the 
CE5 impact glass beads 

A total of 32 CE5 impact glass beads with smooth exposed surface 
(Fig. 1 and Supplementary Fig. 2) and a mare origin (Extended Data 
Fig. 2) were selected for the NanoSIMS measurements of water 
abundances and hydrogen isotopic compositions. In the first 
analytical session, the objective was to investigate the homogeneity 
of water abundance within each glass bead; therefore, two analyses 
of the cores and the rims on relatively large glass beads were carried 
out to identify the hydration and or de-gassing candidate profiles. In 
the second session, six more-detailed profile analyses were carried 
out on five CE5 impact glass beads (Fig. 1), as well as several duplicate 
measurements on other glass beads (Supplementary Fig. 2). Overall, 
the analytical water abundances in the first session are in the range 
of the second session due to the heterogeneous distribution of water 
in the CES impact glass beads (Fig. 2 and Supplementary Table 4). 
The water abundances of the analysed CES impact glass beads vary 
from -0 to 1,909 ug g” and are negatively correlated with their 6D 
values (-990 + 6%o to +522 + 440%) (Fig. 2 and Supplementary Table 
4). Water abundances measured in all CES impact glass beads display 
higher values at the rims and decrease towards the glass bead cores 
(Fig. 3 and Extended Data Figs. 4 and 5). The 5D values display the 
reverse trends (Fig. 3 and Extended Data Figs. 4 and 5). Three profile 
measurements (Profile 3, Profile 4 and Profile 5) show a sharp decrease 
in water abundance from the rim towards the core (Fig. 3 and Extended 
Data Figs. 4 and 5), while the other three (Profile 1, Profile 2 and Profile 
6) display a modest decrease of water abundance towards the cores 
(Fig. 3 and Extended Data Figs. 4 and 5). It is also noticeable that the 
latter profiles have higher water abundances (1,055-1,909 ug g”) 
at the rim than do the former profiles (176-1,423 ug g”) (Fig. 3, 
Extended Data Figs. 4 and 5 and Supplementary Table 4). All of the 
profile measurements acquired on the central areas of CE5 impact 
glass beads have water abundances comparable to the instrument 
H,O background (Supplementary Tables 3 and 4). 


Data availability 

All geochemical data generated in this study are included in 
Supplementary Tables 1-5 and are available via Zenodo (https://doi.org/ 
10.5281/zenodo0.7660603). 
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Extended Data Fig. 1| MgO/AI,0; versus Ca0/Al,0,; for the studied CE5 are from MoonDB (http://search.moondb.org) and are shown for comparison. 
impact glass beads. Most smooth and homogeneous CE5 glass beads plot in Compositions of the CES bulk basalts and bulk soils (shaded yellow square) are 
the mare domain (Supplementary Table 2). Literature data for lunar impact and from Tian et al.” and Li et al.**. Classification of lunar glasses (dashed lines) is 


pyroclastic glasses in Apollo samples, Luna samples, and lunar meteorites (LMs) from Refs. ’”°°. 
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Extended Data Fig. 2 | The water abundance calibration lines established analytical uncertainty is 1.7% (20) for the 1* session and 2.0% (20) for the 2"¢ 
on standards. The linear slopes of the two analytical sessions are comparable session. The error bars are 20. KOV: Kovdor apatite, DAP: Durango apatite, SCOL: 
and consistent with the slope of the instrument’s long-term stability”. The San Carlos olivine. 
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Extended Data Fig. 3 | Hydrogen isotope analysis sequences for the 

two analytical sessions. Analysis sequences are shown to monitor the 
reproducibility of hydrogen isotope analysis of the standards. 6D values of 
Kovdor apatite (KOV), Durango apatite (DAP), SWIFT MORB (MORB), 1833-11, 
and 519-4-1 have been corrected for instrumental mass fractionation. And 

the fractionation factor (Qm) was established by SWIFT MORB glass. The 
reproducibility of 5D analysis throughout the whole analytical session-1 and -2 
are 58 %o (N = 38) and 117 %o (N = 11) (20), estimated by repeated analyses of the 
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SWIFT MORB glass standard. The average SD values calculated for the Kovdor 
and Durango apatite are -39 + 37 %o and -80 + 134 %o (Supplementary Table 3), 
respectively, consistent with their recommended values within analytical 
errors™%”. The average SD values measured for the basaltic glass standards 
1833-11 and 519-4-1 are -54 + 42 %o and -80 + 28 %o (Supplementary Table 3), 
respectively, which is consistent with the 6D values reported in our previous 
work®. The error bars are 2 SD. 
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Extended Data Fig. 4 | Two zoning profiles of water abundances and 6D 
values for CE5 impact glass beads (profile-3 on CE5#33,003 and profile-4 
on CE5#33,046). a. water abundance profile on CE5#33,003. b. 6D profile 
on CE5#33,003. c. water abundance profile on CE5#33,046. d. 5D profile on 
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CE5#33,046. The analytical locations of profile-3 on CE5#33,003 and profile-4 
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on CE5#33,046 are shown in Fig. 1. Hydroxyl/water diffusion modelling results 
at 360 K are shown in black lines. The diffusional durations vary from1to 4.5 
years an 2 to 6 years for profile-3 and —4, respectively. For modelling details see 


Methods. The error bars are 2 SD. 
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Extended Data Fig. 5 | Another two zoning profiles of water abundances 
and ôD values for two CE5 impact glass beads (profile-5 on CE5#33,052 and 


n CE5#33,076). a. water abundance profile on CE5#33,052. b. 8D 


profile on CE5#33,052. c. water abundance profile on CE5#33,076. d. 6D profile 
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on CE5#33,076. The analytical locations of profile-3 and profile-4 are shown in 
Fig. 1. The diffusional durations vary from 0.5 to 4 years and 2 to 14 years at 360 K 


for profile-5 and —6, respectively. For modelling details see Methods. The error 
bars are 2 SD. 
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Extended Data Fig. 6 | Raman spectra. a,b, MORB glasses MRN-G1 (a) and glass bead (CE5#33, 076) measured by NanoSIMS was chosen for Raman profile 
EPR-G3 (b). c, CE5 impact glass bead, CE5#33,076. Raman peak at -3500 cm ‘is analysis. The Raman analytical positions of CE5 impact glass beads are outlined in 
ascribed to hydroxyl/molecular water. The grey lines are apparent baseline of Supplementary Fig. 2. 


the Raman peak between -3300 and 3700 cm”. A notable H,O zoning CE5 impact 
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